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Introduction 



A methane-based meteorological cycle involving clouds and rain has been suspected 



to occur on Titan since the time of Voyager, when methane was measured by IRIS (Hanel 



et al. 1981). Pioneering radiative transfer models of the atmosphere (Toon et al. 1989 



McKay et al. 1989) were used to interpret these spectra, predicting rain and patchy clouds 



from 10 - 30km (Toon et al. 1988). Ground-based measurements at shorter wavelengths 



supported the presence of clouds and discovered their daily variation (Griffith et al. 1998 



2000). Advances in imaging technologies lead to spatially-resolved observations of clouds 



near the south pole (Brown et al. 2002 Roe et al. 2002), as well as a region of clouds 



confined to southern mid-latitudes (Roe et al. 2005a). Since the arrival of Cassini it has 



been possible to resolve complex cloud structures (Porco et al. 2005), measure updraft 



velocities (Griffith et al. 2005), and identify polar clouds that may be composed of ethane 



(Griffith et al. 2006). 



Global scale precipitation has been suggested using the methane relative humidity 
profile measured by Huygens ( Tokano et al.||200"6 ). While it was first predicted that droplets 



would evaporate before reaching the surface (Lorenz 1993), recent models suggest rain 



can fall to the ground (Graves et al. 2008). Fluvial channels indicate that a liquid was in 



contact with the surface at some time (Tomasko et al. 2005). 



Near-IR observations in two methane transmission windows (1.6 and 2.0 /im) have been 
used with radiative transfer models to measure the column of condensed-phase methane on 



Titan (Adamkovics et al. 2007). These data revealed a global cloud of solid methane in the 



25 - 35 km altitude range, with an enhancement of opacity on the morning hemisphere that 
is consistent with methane drizzle. Recently, these results have been called into question 



using observations made at the Gemini-North telescope with the NIFS spectrometer (Kim 
et al. 1120081). 
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Here we present measurements over consecutive nights, together with radiative transfer 
modeling, to show that the signature of increased condensed-phase methane cannot be 



due to surface reflectivity artifacts as suggested by Kim et al. (2008), and use high 
(Sun-target-observer) phase angle measurements to rule out phase angle artifacts. We 



account for the conclusions and analysis of Kim et al. (2008) by quantifying the errors 
that can arise in the surface subtraction due to noise and bias. Our observations indicate 
that while precipitation occurs predominantly near Xanadu, it is not solely a morning 



phenomenon as suggested earlier (Adamkovics et al. 2007). 



2. Observations 

Ground-based observations were performed between 2006 Dec 28 and 2007 Mar 11 



UT with the adaptive-optics aided integral-field spectrometer, SINFONI (Eisenhauer et al. 



2003), at the Very Large Telescope (VLT). The spectrometer uses two sets of 32 stacked 
mirrors to divide the field of view (FOV) into 32 'slitlets', which are then optically arranged 
into a single synthetic long slit that is dispersed onto a cryogenically-cooled 2K x 2K 
Hawaii detector. Each of the 32 slitlets is imaged onto 64 x 2000 pixels on the detector. 
We used the 0.8" x 0.8" FOV, corresponding to a pixel scale of 0.0125" x 0.025", with the 
H+K grating covering 1.45 - 2.45 jum at a resolving power, R = A/AA, of 2000 to 3400. 



Data-reduction follows the SINFONI pipeline (Modigliani et al. 2007). The pipeline 
includes a reference bad-pixel map, a check for pixels approaching a non-linear detector 
response, dark current correction, and flat fielding. All of the spectra are distortion 
(curvature) corrected using standard flats and arc lamps frames, along with wavelength 
calibration from a reference arc line table. Spectra are interpolated onto a common 
wavelength grid before reconstruction of the datacube using the slitlet position and slitlet 
edge tables. 
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Mosaicking is performed using AO pointing keywords to overlap images to the nearest 
pixel and then taking the median value for overlapping pixels. Telluric correction is 
performed with the extracted spectra of standard stars. Flux calibration is performed by 
normalization to the disk-integrated stratospheric flux from Keck/OSIRIS observations 



(Adamkovics et al. 2007), and are reported in units of albedo, I/F. The observation dates 
and viewing geometries are given in Table 1. 

The Visual and Infrared Mapping Spectrometer (VIMS) on the Cassini spacecraft uses 



a spectral mapping technique (Brown et al. 2004) to assemble images at 352 wavelengths, 



covering 0.3 to 5.1 /im. Channels near 1.6 and 2.0 /im can be used for direct comparison 
with ground-based observations, and provide the opportunity for high phase angle viewing 
geometries that are not obtainable from Earth. We use VIMS data that are publicly 
available from the T10 flyby on 2006 January 15 UT. The data are reduced from observed 



data numbers to units of albedo, I/F, according to Barnes et al. (2007) 



Radiative Transfer Model 



In order to interpret the observed spectral datacubes, we create a forward model 
to simulate the flux that is reflected from Titan's atmosphere and surface. Our model 
incorporates well-established numerical solutions to the radiative transfer (RT) equation 



(Griffith et al. 



1991 



Brown et al. 



2002 



Adamkovics et al. 



2004 



2007) that reproduce 



near-IR observations exceptionally well. Updates to the model presented in Adamkovics 



et al. (2007) include the treatment of the surface reflectivity; including the specification of a 



unique surface spectrum for Xanadu, averaging the surface albedo map over the footprint of 
a pixel in the model, and an explicit fit for converting the observed surface albedo map to a 
surface reflectivity map that is input to the model (described in Section 3.3). Limitations 
in current analyses are due to input parameters (e.g., CH 4 opacities, surface albedo spectra, 
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and aerosol scattering) and not the numerical methodology or assumed geometry. A 
2-stream approximation (at small phase angles) or discrete ordinates method (at high 
phase angles) is used to solve the radiative transfer equation for 16 pseudo-plane-parallel 



layers from - 200 km altitude (Toon et al. 1989 McKay et al. 1989 Stamnes et al. 



1988). We correct for the curvature of the atmosphere with an established geometrical 



correction (Tran and Rannou 2004). We use the atmospheric temperature and pressure 
profiles measured in situ by Huygens/HASI ( |Fulchignoni et aL 2005), which are in excellent 



agreement (below 200 km) with the Voyager profiles (Lellouch et al. 1989). Starting with 



the in situ measurements of aerosol extinctions made by Hugyens/DISR (Tomasko et al. 



2005), we vary the extinction at other locations on Titan to reproduce our observations of 



the global distributions of haze. The uncertainty of the CH 4 opacity, which is greatest for 
weak transitions, translates into an uncertainty in both the retrieved surface albedo and to 
a lesser extent the tropospheric aerosol extinction that is constrained at other wavelengths. 



3.1. Gas opacity 



The vertical profile of CH 4 mixing ratio measured by Huygens/GCMS (Niemann 



et al. 2005) is applied globally and the absorption coefficients for CH 4 presented by Irwin 



et al. (2005) are used to calculate CH4 opacities. Corrections to these coefficients have 



been published using the measurements made by Huygens/DISR (Tomasko et al. 2008c), 



however the corrections do not extend to 2 /im and are not applied here. An alternative 
determination of the CH4 opacity may be performed with line-by-line calculations that use 



the very high resolution, high sensitivity, CH 4 spectra of Brown (Brown 2005). These data 



are currently available in the HITRAN database. However, the line strengths for these 
transitions are only applicable at room temperature, and quantum mechanical assignments 
of individual CH 4 lines are required before opacities can be calculated at the temperatures 
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applicable to Titan. As described by Brown (2005), assigning these transitions requires 



considerable laboratory and computational effort. 

Other sources of gas-phase opacity in the 2.0 /zm band include low temperature 



collision- induced absorption by N 2 and H 2 (McKellar 1989). The H 2 -N 2 collision complex 



is a significant source of opacity near 2.1 /im, which decreases the observed slope of Titan's 
spectrum in the 2.0 /im band relative to the 1.5 /im band. Since the magnitude of the 
observed flux can be roughly mapped onto the altitude where most of the scattered photons 
originate — that is, regions of weak atmospheric absorption are bright and probe the 
surface, whereas regions of high gas opacity are dim and probe the upper atmosphere — the 
gradual variation in the 2.0 /xm band extinction facilitates altitude discrimination relative 
to the 1.5 /im band. This collision-induced opacity is critically dependent on the mixing 
ratio of molecular hydrogen, which we take to be well mixed throughout the atmosphere at 



levels of ~0.1% (Samuelson et al. 1997). 



3.2. Aerosol Haze Scattering 

Observations at IR wavelengths can be represented by models with an approximated 
scattering phase function (e.g., a one-term Henyey-Greenstein). When calculating the IR 
flux reflected from Titan at low phase angles there is, in general, a degeneracy between the 
aerosol scattering phase function and the extinction. For example, an increase in albedo 
can result either from a greater aerosol density or aerosols that are more backscattering. 

The DISR instrument on Huygens was used to measure the haze extinction profile 



from 0.5 - 1.5 fim (Tomasko et al. 2005), and our ground-based observations made during 
probe entry (at 1.5 /mi) use these measurements as a benchmark for the extinction near the 
latitude of the Huygens landing site, 10°S. The cumulative DISR extinction is extrapolated 
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to observations at 2 /iin (Adamkovics et al. 2007). Recently, a parameterization of aerosol 



scattering in three altitude regimes on Titan has been fit to the DISR data (Tomasko 



et al. 2008a). The authors propose a model of the wavelength dependence of the scattering 



phase function that is derived using the DISR data and scattering that is treated with the 



T-matrix method (Mishchenko et al. 1996). The net solar flux both in Titan's atmosphere 



and at the surface, calculated using the detailed optical properties and distribution of 



particles, is in good agreement with the values calculated by McKay et al. (1989) (Tomasko 



et al. 2008b). 



We use a simplified parameterization of the phase function, following McKay et al. 



(1989). The two-stream methodology we use incorporates the total flux that is scattered 



in the forward and backward directions. This method is not sensitive to further subtleties 
in the scattering phase function. For the 1.5 and 2.0 /zm windows, a one-parameter 
Henyey-Greenstein phase function has been adequate for interpreting the observations and 



has accurately predicted the heat balance in the atmosphere (Tomasko et al. 2008b) 



The spatial variation in aerosol density is modeled using the method presented in 



Adamkovics et al. 



(2006). The stratospheric aerosol extinction increases northward at a 



rate of 0.55% per degree latitude from 45°S up to ~70°N, the highest latitudes observed 
here. The tropospheric aerosol extinction is uniform below 40 km altitude. 



3.3. Surface reflectivity 



The surface reflectivity is modeled with the simplest representation that can reproduce 
both spatial and spectral variations in the observations. We start with a high-resolution 



cylindrically-projected map of the 2.018/mi albedo (Barnes et al. 2007) and re-project it 



onto a sphere viewed on the plane of the sky, at the instrument plate scale, using the 
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JPL HORIZONS ephemerides. For each detector pixel the albedo map is integrated over 
the spatial region that corresponds to the polygon enclosed by the corners of that pixel, 
creating an albedo array, a x>y . If the corner of a pixel is located off the edge of the disk 
then the region of the albedo map that is enclosed by the remaining three corners is used. 

Since the observed 2.018 fim albedo has a contribution from atmospheric scattering, it 
would be a mistake to use a X)V as the input for the surface reflectivity, R x , y , in the model. 
This would result in an overestimate of the flux at surface-probing wavelengths. We use a 
linear scaling and offset of the entire map to convert from albedo to reflectivity, 



Rx,y = Oix,yPl + Vl (!) 

where the parameters are determined by fitting the output model spectrum to the 
observations. A series of models are calculated over a range of input values for pi and p2- 
In each case the mean squared residual between the model and observations is calculated 
over the entire datacube, and the parameter combination corresponding to the model with 
the smallest residual is selected as the best fit, pi=0.83 and p2=0.02 

In the radiative transfer calculation we assume isotropic scattering at the surface. As 



described in Toon et al. (1989), the total upward flux from the surface is the sum of the 
reflected downward diffuse flux and flux from any other sources, S B fc- At 2 fim, the other 
relevant source of flux is the reflection of the incident sunlight that reaches the surface, 
given by 

SWc = Rx, y ^o exp {-t/hq)'kF s (2) 

where /x is the cosine of the solar incidence angle, r is the optical depth of the atmosphere, 
and 7rF s is the incident solar flux. 

The Huygens probe measured the surface reflectivity spectrum in situ from the 
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visible up to 1.5 /im (Tomasko et al. 2005). At longer wavelengths, remote observations 
through narrow regions of low methane opacity are only sensitive to portions of the surface 
reflectivity spectrum. Since a flat spectrum does not reproduce the observations, i.e. the 
surface reflectivity in the 1.6 /im window is systematically larger than in the 2.0 /im window 



and the contrast between bright and dark regions is different ( Coustenis et al. 2005 de 



Pater et al. 2006), we incorporate a surface reflectivity spectrum in our model. 



We use a set of Gaussian features in the surface spectrum to resolve discrepancies 
between the calculated albedo and the observations. The surface reflectivity for the model 
datacube is given by 

n 

R(x,y,X) = J R,,,J](l-A l exp(-(A Cii -A)/5A i ) 2 ) (3) 

where we call the product term the surface reflectivity scale factor. Each spectral feature, 
i, is defined by a central wavelength, A Cji , line width, <5Aj, and amplitude, A iy for n features. 
Values for the Gaussian features are given in Table [2} Three narrow absorptions are used 
in the base model and (since R x ,y is determined at 2 /im) another Gaussian feature is used 
to represent the increase in reflectivity toward 1.6 /xm. It is unclear if these transitions can 
be uniquely related to the surface composition. The surface reflectivity scale factor is fit 
globally and individual spectra in the observed datacubes are integrated over vast regions 
covering more than 10 5 km 2 . 

The surface reflectivity scale factor is plotted in Figure [l] with an example comparison 
of the model output and observations. By considering the surface reflectivity a free 
parameter, it is possible to mask uncertainties in atmospheric opacity and fit the spectra 
arbitrarily well at wavelengths that sense the surface. The relative contribution of surface 
flux depends on wavelength, surface reflectivity, and viewing geometry. To approximate 
the fraction of flux arising from the surface in the 2 /im window, we compare a spectrum 
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calculated with the base model to a model spectrum with the surface reflectivity, R x>y , 
set to zero; the two spectra are plotted in Figure |2j The minimum contribution of the 
atmospheric flux is 17% of the total , at 2.04 /zm for a surface reflectivity of 0.12, and will 
be larger for lower reflectivity surfaces, and locations near the edge of the disk where the 
surface is limb-darkened while the atmosphere limb brightens. 



3.4. Condensed-phase methane 

After Voyager it was suggested that methane is saturated in some regions of Titan's 



'Flasar et al. 


1981; 


Toon et al. 


1988 



measured by Huygens (Fulchignoni et al. 2005) and the presence of clouds both conclusively 



indicated that methane exists in solid and liquid form in Titan's atmosphere. Despite this 
evidence, RT models used for direct comparison to near-IR observations have only recently 



included a rudimentary treatment of condensed-phase methane opacity (Adamkovics et al 



2007). The morphology of condensed-phase methane (e.g., particles, droplets, or grains) is 



unconstrained by these models. A solid cloud of methane (Tokano et al. 



2006 



Adamkovics 



et al. 2007) is included in the 25-35 km altitude layer of the model with a methane column 
volume of 1.5 cm 3 cm" 2 . The column volume is equivalent to a uniform layer of condensed 
methane with a thickness, £=1.5 cm. The optical depth of condensed-phase methane is 
given by 



r c (A) = £a(X) 



(4) 



where the absorption coefficient spectrum, a(X), is from laboratory measurements (Grundy 
et~aLl[2002| . The base model has no additional methane precipitation. This basic treatment 
of the condensed phase opacity assumes only absorption and negligible scattering — i.e, the 
scattering at these altitudes is dominated by aerosol haze. 



3.5. Point spread function 



The surface reflectivity (from the VIMS albedo map) and the stratospheric aerosol 
gradient (a free parameter) are the two spatially-varying inputs in the model. After these 
inputs have been determined, the RT calculation of a spectrum is performed for each of 
the approximately 4000 spatial locations, building a model datacube at the plate scale 
of the observations. The last step is to convolve the model with a point spread function 
(PSF). We use the telluric calibration star datacube, collapsed in wavelength across H- and 
K-bands separately, as the reference images for the PSF. A 2d Gaussian is fit to the PSFs 
and used as the convolution kernel for each band. Fitting the PSF and then using centered 
convolution kernel avoids the addition of noise into the model. 



4. Surface Subtraction 

Atmospheric scattering at 2 /mi accounts for at least ~17% of the total observed flux, 
the rest is from the surface. Changes in surface reflectivity at different locations on Titan 
generally swamp the signal from spatial variations in atmospheric scattering. 

Images probing the lower atmosphere are also sensitive to the surface. The latter can 



be subtracted as described below, to discern atmospheric phenomena (Adamkovics et al. 



2007). We calculate surface subtracted, AI/F, images using the method of 



Adamkovics 



et al. (2007) for all of the observations and corresponding models. Each pixel in the 



subtracted images is defined according to 



AI/F = (I/F) Xl -fx(I/F) M , (5) 

where the bandpasses for (I/F)\ 1 and (I/F)\ 2 are given in Table [3] The bandpass for 
(I/F)x 2 is selected to cover the wavelength region that probes the surface with the least 
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atmospheric contribution, while (I/F)\ 1 is integrated over a bandpass with increased 
methane absorption and therefore contains a larger contribution from the atmosphere. 
The linear Pearson correlation coefficient, r, between the AI/F image and the (I/F)\ 2 
is determined for a range of /. Surface subtracted images using a range of values for / 
demonstrate that small values (/ < 0.57) yield images correlated with surface probing 
images, while large values (/ > 0.74) are anti-correlated with the surface image — that is, 
they yield and inverted image of the surface albedo pattern and large negative r (Figure [3]). 
Intermediate values minimize the absolute value of the correlation coefficient. Although 
the contribution of flux from the surface depends on a reflectivity, the value of / that 
corresponds to minimum |r| is the optimal value used in this simplistic image subtraction. 
A single / may be used in a model-independent, approximate method for removing surface 
albedo variation from images. 

Since the background noise in regions of the image with sky can affect the statistical 
analysis, we performed the calculation of / for each night using only the pixels on the disk 
of Titan, yielding /=0.67, with a mean absolute deviation cr=0.011. The scale factors 
for each night are presented in Table [4] for the 1.5 /zm and 2 /im bands. These values are 
significantly lower than when the analysis is performed on the entire image, which includes 
noise from the sky; in that case /= 0.76 and cr=0.013. The values of / calculated while 



including regions of the sky are consistent with the smallest value calculated by Kim et al 



(2008) and larger than /=0.72 used in Adamkovics et al. (2007). The AI/F images for the 



observations and all models — each calculated with /=0.67 — are presented in the bottom 
four rows of Figure |4} 

Comparison of the observations and the models, in both wavelength regions (Ai and A2) 
demonstrates that where the surface albedo is well mapped by Cassini/VIMS, the models 
reproduce the observations well. Additionally, changes due to additional condensed-phase 
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methane opacity do not produce significant changes to the images in the A2 bandpass. 
The small variations in the flux from subtle changes in opacity and reflectivity are more 
clearly observed after subtraction of the surface albedo variation. The la pixel-to-pixel 
noise is between 0.0012 and 0.0016 for all of the observed AI/F images. These values are 
calculated two ways: first, a conservative estimate is made by calculating the mean of the 
absolute value of the residuals between the base model and the observations, and secondly 
by calculating the FWHM of the histogram of the AI/F values at the center of the disk in 



each image. In the following discussion, we investigate the suggestion of Kim et al. (2008) 
that bright regions of the surface cause artifacts in the difference images. 



5. Discussion 

The ground-based observations and RT models are presented in Figure |4j Although the 
entire observed datacube is calculated (covering 1.45 - 2.45 /mi) we focus here on images 
from the wavelength windows that are used to identify condensed phase methane opacity. 



Kim et al. (2008) argue that anti-correlations are created in AI/F images due to 
regions of high surface albedo having different / values than dark regions. The argument 
predicts that bright regions will, in general, be anti-correlated with dark regions in AI/F 
images. Observations in Figure [4] demonstrate that a gross anti-correlation does not occur. 
The brightest regions in the top row of images does not correspond with the dark regions 
in the subtracted images. Due to limb darkening, regions of the surface with the brightest 
(I/F)\ 2 always occur at the center of the disk, whereas dark regions in the AI/F (if present 
at all) occur near the limb. There is no gross anti-correlation between (I/F)\ 2 and AI/F 
images. 



Perhaps the argument of Kim et al. (2008) applies only to the region of highest surface 
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reflectivity, Xanadu, or as they suggest, the spectrum of the surface needs to be taken into 
consideration. The four consecutive nights of observations from 2007-01-28UT through 
2007-01-31UT are used to test this hypothesis. On 2007-01-28UT the eastern tip of Xanadu 
is visible on the limb of Titan, and by 2007-01-31UT this location is observed at the center 
of the disk. In the AI/F images from these dates, there is a darker region over the tip of 
Xanadu on 2007-01-28UT, however, the same location is no longer dark on 2007-01-31UT. 
This sequence of observations suggest that artifacts related to static properties of the 
surface — which are expected to track with rotation — do not account for the dark regions 
of the surface. To further test this hypothesis, we employed a systematic change in the 
surface reflectivity spectrum, but only near Xanadu. 

The Xanadu region is outlined (somewhat arbitrarily) using a high resolution ISS 



map (Porco et al. 2005). A test case is modeled where the spatial locations in the 
observed datacubes that are bounded by Xanadu (plotted in the top row of Figure [ij 
have two additional surface absorption features, listed in Table |2j and shown graphically 
in Figure [TJ The features are centered on the wavelength regions that are used to probe 
the atmosphere, above the surface in both H- and K-bands, and are narrow enough to fit 
within the A2 bandpass used for the images. The effect of this mildly contrived change 
on the surface-subtracted images is demonstrated in the bottom row of Figure |4j Not 
unexpectedly, the region of Xanadu is darker in the difference images. Altering the surface 
spectrum in such a way could be modified slightly to fit an individual observation very 
well. However, the series of observations from 2007-01-28 through 2007-01-31 illustrate 
that a systematic change to the surface reflectivity over Xanadu cannot reproduce all the 
observations. For these dates, the mean of the absolute value of residuals between the 
condensed-phase model and the observations is smaller than the residuals of the model with 
the artificial surface spectrum. On 2007-02-23, 2007-03-10, and 2007-03-11, the models 
with either additional methane opacity or the surface spectrum change both have larger 
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residuals than the base model. 

The observed AI/F image on 2007-01-29 has a darker region near the limb than in 
the model with a unique Xanadu surface spectrum. This suggests that the absorption 
spectrum of the Eastern tip of Xanadu is underestimated in the model. However, on 
2007-01-31 this same region of Xanadu is near the center of the disk and the absorption 
spectrum is overestimated in the model. A change to the surface reflectivity spectrum 
cannot simultaneously match both observations. A similar effect was noted qualitatively in 
(2007), and is also illustrated in Figure [i] with the observations from 



Adamkovics et al. 



Kim 



et al. (2008). Whereas the two regions of Eastern and Western Xanadu are both labelled £ A' 



in Figures 1 and 3 by Kim et al. (2008), here we label the East and West regions separately 
as 'A' and 'Al', respectively. This clarification shows that while 'Al' is a bright albedo 
region that shows up as a dark feature in the Keck/OSIRIS difference images from 2006 
Apr 17, this same bright feature — near the center of disk — is no longer a dark region in 
the difference image viewed with VLT/SINFONI. 

We test for surface scattering phase artifacts by performing our analysis on 
Cassini/VIMS data from the T10 fryby. Since the dark regions in the AI/F images occur 
near the limb, it is possible that some property of the surface scattering is the cause of 
an artifact that causes Xanadu to always be dark in the AI/F images near the limb. 
However, the images in Figure [6] demonstrate that Xanadu is not an inherently dark region 
in AI/F images. Intensity variation in the VIMS AI/F image does indicate that a single 
value of / will not remove all the spatial variation in the subtracted image, however, the 
gross over-correction for bright regions is avoided for a minimum value of the correlation 
coefficient, |r|. Our RT models also reproduce the high phase angle observations with 
Lambertian surface scattering and a surface reflectivity scale factor spectrum that does not 
vary across the disk. 
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The observations and models presented here can be reconciled with the measurements 



of Kim et al. (2008) by quantitatively evaluating factors that can affect the determination 
of /. We calculate / for a subset of pixels in a circle centered on Titan and plot the 
resultant / as a function of Titan radius, Rt, in Figure [7j. Near the center of the disk 
(Rt < 0.7) , / is roughly constant. There is a decrease in / when limb-darkened regions - 
where contributions from the surface are small — are included in the correlation analysis. 
Including regions of the sky adds noise and increases the calculated value of /. We further 
investigate this effect by systematically adding noise to the images and reproducing our 
analysis. We generate images with Gaussian noise (of a given FWHM in I/F), add the 
image of noise to our observations, and then recalculate / (lower scale bar in Figure [7]). 
Additional noise leads to a systematically higher determination of / (dotted line Figure [7J) 
such that 1% of additional noise can lead to a increase in / by ~0.02. A linear bias 
is similarly tested by scaling the (I/F)\ 2 image (dashed line in Figure [7]). This is an 
approximation of uncertainties due to flat-fielding, bias-correction, or telluric correction. 
Both positive or negative deviations in / can occur due to biased images. Data from 
2007-01-28UT are used above for illustration, however the results generalize to each night 
of our observations. These tests demonstrate that using the entire image in the analysis, 
an offset during data reduction, or a noisier dataset can all lead to larger values of /, and 
thus an 'over-subtraction' of the surface. Overestimating / will lead to anti-correlated 
artifacts from regions of bright surface reflectivity, as detailed in 



Adamkovics et al. 



(2007), 



and evidenced in Kim et al. (2008). Figure M demonstrates that the Gemini/NIFS data 



presented by Kim et al. (2008) are noisier than either the VLT/SINFONI or Keck/OSIRIS 
images presented in Adamkovics et al. (2007). 



-17- 



6. Conclusions 

Atmospheric phenomena are most likely responsible for the observed dark regions in 
the AI/F images and have been quantitatively reproduced in models using additional 



condensed-phase methane (as in Adamkovics et al. (2007)) and tested here against 



surface scattering artifacts. Since the extinction properties of haze are not strongly 
wavelength-dependent in the 2 jum region, spatial variations in haze are removed during the 
image subtraction and cannot reproduce the dark regions in AI/F images. Altering the 
haze properties enough to create the dark AI/F regions results in discrepancies that fail to 



reproduce the observed albedo in tropospheric or stratospheric images (Adamkovics et al. 



2007). Localized meteorology has already been observed (Roe et al. 2005b), so perhaps 
precipitation is also geographically linked. On 2007-03-09 we observe condensed-phase 
methane over Xanadu on the evening limb, Figure [4j indicating that precipitation can occur 
either in the morning or in the afternoon. 

Further modeling of scattering by strongly absorbing droplets or ices is needed to 
determine if the rain can reach the ground and if there is a predicted signature of rain at 
longer wavelengths (e.g, 5/im). There is also a possibility that the condensed methane 
coexists as a surface coating on aerosols. Alternatively, the optical properties of haze may 
be such that they mimic methane transmission spectrum, which is unexpected but cannot 
be ruled out. Currently, enhancement in condensed-phase methane over Xanadu remains 
the least speculative and the only modeled description of dark regions in AI/F images. 
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Table 1: Observation dates and viewing geometry. 

Obs. Date Airmass Diameter Sub-observer Point Sub-solar Point Phase Angle 
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* Sub-solar longitude is greater than sub-observer longitude. For these observations, the 
sub-solar point is to the left of the center of the disk in the images in Figure [1] 



Table 2: Surface reflectivity spectrum features 



Line center, Line width, Amplitude, Notes 
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* Since R x , y in Equation 3 corresponds to observations at 2 /im, 
a negative amplitude is used to represent the increase in reflectivity 
toward 1.6 /xm (Figure [lj. 
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Table 3: Bandpasses (jum) used for AI / F images 
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"Approximated bandpass for individual VIMS channels centered at 
2.034, 2.068, 2.117, and 2.150 /zm. 



Table 4: Surface subtraction parameters. 

Obs. Date (UT) 1.5 /xm band 2.0 /im band 
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Fig. 1. — Surface reflectivity spectrum used in RT models is a combination of the surface 
albedo measured by Cassini/ VIMS at 2.018 /im and a wavelength-dependent scale factor 
(dotted line, left axis). The surface reflectivity spectrum is used to calculate the I/F spec- 
trum (thick line, right axis) and compare to observations (thin line, right axis). The base 
model surface reflectivity scale factor (dotted line, left axis) fits both the dark and bright 
regions, and a contrived spectrum with two additional absorptions (indicated near '*') is 
used only for the region of Xanadu to test the differential imaging analysis. 
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Fig. 2. — The contribution of flux from the surface in the 2 /xm window was calculated by 
comparing a spectrum of a bright region near the center of the disk in the base model (solid 
line) with a model that has the surface reflectivity, R x , y , set to zero (dotted line). In the 
case that R x ^ y =0, all of the flux is from scattering in the atmosphere. The atmospheric 
scattering is assumed to be the same in the base model, which also includes the light that 
is reflected from the surface. The maximum surface contribution occurs at 2.04 /mi, where 
the flux from the atmosphere is 17% of the total flux. Integrating the surface contribution 
over the entire 2 /im window gives a minimum atmospheric contribution of 32% of the total 
observed albedo. 
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Fig. 3. — The scaling factor for the image of the surface flux is determined by an empirical 
test of the correlation between the surface probing image and the subtracted image. 
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Fig. 4. — K-band (2 jttm) VLT/SINFONI observations and radiative transfer models both 
with and without an increase in condensed phase methane near Xanadu. The increase is 
modeled to occur from 30°S to 20°N, and from the limb to within 35° of the sub-observer 
longitude for all observations. A prominent southern mid-latitude cloud (identified at longer 
wavelengths) contaminates the AI/F image on 2007-01-30 UT, appearing as a bright feature 
between 30 - 45° S. An incorrectly high value for / affects areas with the same surface albedo 
uniformly and would result in all of Titan's bright regions showing small AI/F values. This 
is not the case here. While Xanadu corresponds to low AI/F, other high-albedo surface 
terrains such as Dilmun (10°N 175°W), Tsegihi (40°S, 35°W), and Adiri (10°S 210°W) do 



not show similar negative excursions — however, this is not the case for Figure 3b in Kim 



et al. (2008). 



Fig. 5. — Comparison of field- integral spectrometer images from three instruments: 
VLT/SINFONI, Keck/OSIRIS, and Gemini/NIFS. The VLT and Keck data are from 



Adamkovics et al. 



(2007), while with lower S/N and lower resolution Gemini data are com- 



bined from Figures 1 and 3 in Kim et al. (2008). For clarity, the two regions labeled in 'A' 



by Kim et al. (2008) are separated here as 'Al' and 'A'. The 'Al' regions are dark in the 
difference image from Keck while appearing bright in the VLT observations. 




Fig. 6. — Cassini/VIMS observations and surface subtracted images at high phase angle 
(top), with radiative transfer models of the same observing geometry (bottom). The outline 
indicates Xanadu. Limb brightening in the AI/ F images due aerosol scattering is reproduced 
in the models. Bright scattering at the poles is due to optically thin clouds. 
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Fig. 7. — The optimal surface-subtraction scale factor / depends on the region of interest in 
the image (solid line, lower axis), bias between the (I/F)\ 1 and (I/F)\ 2 image (dashed line, 
top axis), and the noise in the analyzed images (dotted line, lowest axis). Calculations for 
bias and noise were made using Rt < 0.7. 



